Exactly Solvable Ginzburg-Landau theories of Superconducting Order Parameters 

coupled to Elastic Modes 
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We consider two families of exactly solvable models describ- 
ing thermal fluctuations in two-dimensional superconductors 
coupled to phonons living in an insulating layer, and study the 
stability of the superconducting state with respect to vortices. 
The two families are characterized by one or two supercon- 
ducting planes. The results suggest that the effective criti- 
cal temperature increases with the thickness of the insulating 
layer. Also the presence of the additional superconducting 
layer has the same effect. 

PACS numbers: 05.30-d, 74.20De, 11.10-z 



I. INTRODUCTION 



Recently Fateev introduced two families of in- 

tegrable models which can be interpreted as Ginzburg- 
Landau free energy functionals describing thermal (clas- 
sical) fluctuations in superconducting films. According to 
this interpretation the models of the first family, which 
we shall further call type 1 models, describe a single su- 
perconducting layer deposited on an insulating substrate 
consisting of n layers. The superconducting order pa- 
rameter interacts with the elastic modes of the substrate. 
The effect of the interaction is to shift the local transi- 
tion temperature. Type II models describe the situation 
of a double layer where an insulator is sandwiched be- 
tween two superconducting films. The latter ones may 
have different critical temperatures. 

These models provide a rare opportunity to go beyond 
weak coupling description and obtain non-perturbative 
results for layered superconductors interacting with an 
insulating stratum. In the case of type-I models our 
principal interest is to find how the Kosterlitz-Thouless 
transition temperature depends on the thickness of the 
substrate. Type-II models allow us to study the effects 
of interaction between the superconducting order param- 
eters. 

The paper is organized as follows. In the next Section 
we describe the Fateev's models and explain why they 
can be interpreted as effective Ginzburg-Landau theories 
for layered superconductors. In Section III we study the 
stability of the type I models with respect to the vortices. 
We do a similar analysis for the type II models in Section 
IV and discuss the results in Section V. 



II. EFFECTIVE GINZBURG-LANDAU THEORY 

In this section we describe the Fateev's models and 
construct the effective Ginzburg-Landau free energy. The 
models of the type I describe a complex bosonic field, 
A, interacting with a n-component real scalar field <j} — 
4>n)- There are different models with slightly dif- 
ferent actions for the scalar fields; for our purposes it is 
sufficient to describe only one of them where the classical 
( Euclidean ) action has the following form: 



(1) 



exp( 



1 /Vf2 



2 J2 exp( 
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<j>i + l) + 2exp{<j>n)]} 



This action corresponds to a complex sinh-Gordon model 
coupled to affine Toda chains. The complete list of the 
models of the two families is given in For small n 
form (^ requires certain modifications, in particular for 
71 = it becomes B : 
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The type II models describe two complex bosonic 
fields, Ai_2, interacting with n elastic modes. Again we 
present only one specific model where the elastic part co- 
incides with the one described by Eq.(|l|). The action is 
the following: 



' L s=l,2 



d^Asd^As , 2m2AiAi (3) 



A.A,, 



exp(-0i)-H2Mo2A2A2 exp(- 



I 2exp(-0i) + 2^ exp(0i - (/)j+i) -I- 2exp 
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For n=:0 the above expression has to be modified as: 
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s=l,2 

{AiAi + A2A2 + 2(AiAi)(A2A2)}] 
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As a consequence of the U{1) (type I models) or U (1) x 
U{1) (type II models) symmetry, the models have the 
following conserved charges (here we work in Minkovsky 
space-time): 

where s = 1 for type I and s = 1, 2 for type II models. In 
presence of external chemical potentials the Hamilto- 
nian is modified: H = Hq — h^Qs- One can introduce 
new variables 



with 



A^e 



ith. 



(6) 



to remove the terms with one time derivative. Then we 
obtain the general form of the Ginzburg-Landau free en- 
ergy for our layered superconductors: 



(7) 



where: 



1 + AA 



+ -(9^</>)2 + 2M2 
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and 



s=l,2 



1 + A,A, ' 2 

-f2Mo2Ai Ai exp(-0i) -t- 2A/02A2A2 exp(-0„) 

n-l 

2exp(-(/)i) + 2^exp((/>i - + 2exp((/)„) 



(9) 



Clearly the fields A^ can be interpreted as a supercon- 
ducting order parameters coupled to the optical phonon 
modes 0„. We note that the elastic modes are not 
harmonic, which is necessary to make the model inte- 
grable, but in the limit 7 — > we reproduce the conven- 
tional elcctron-phonon interaction. Indeed, rescaling the 
phonon field (j) — )■ 70, we see that the limit 7 — *■ corre- 
sponds to harmonic phonons. The fact that the spectrum 
of the Toda chain coincides with the spectrum of har- 
monic phonons Q gives us grounds to believe that the 
unharmonicity present in the model does not influence 
the qualitative validity of our results. At the same time 
it is interesting to have a model with phononic spectrum 
that is not restricted to the harmonic one. 

In order to get a better understanding of the nature of 
the superconducting state, let us have a deeper look on 
the free energy for the case n — 0. For the type I models 
the free energy can be written as follows: 



1 + AA 



V{A,A) 



(10) 



V{A,A) 



M^AA 
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For Mq < h the effective potential for the order pa- 
rameter has a minimum at |A| ^0 which signal the ap- 
pearance of superconductivity. Indeed, expanding the 
action around the minimum we notice that the quanti- 
ties T = (Mo2-/l2)/72 

may be interpreted as the distance 
from the mean field transition {T/Tc — 1). Increasing h 
one can go from the disordered to the superconducting 
state. One can see it more explicitly using a semiclassical 
analysis valid for 7^1. Under the transformation 



A = sinhpe*'^ 
the free energy becomes: 



(12) 



^0'^ d^Hid.P? + tanh^ pid,^)' + V,)j{p)} 

(13) 



where: 



Vef f {p, h) = A/q sinh^ p - tanh^ p 



(14) 



The measure of integration in the path integral changes, 
but this is not important for what we want to discuss 
here. For h/2 > Mq, Veff{p,h) has a minimum at 
p — po, where exp2po — ^/2{h/M). Thus when h ex- 
ceeds the threshold a gapless Goldstone mode appears. 
To see it explicitly we expand the free energy around the 
minimum. Rewriting p = po + ^ and keeping only the 
quadratic terms in ^, we get: 



il) 



r 



rf'2;{(a^0' + {Moh/2)e + tanh' Po{d^^f}. (15) 



We can then identify the gapless mode ip, which velocity 
is equal to the bare one, and a massive field ^. 

A similar analysis can be performed on the type-II 
models (again we consider only the simplest case n — 0). 
In this case the transformation 

A« =sinhp«e*'^» (16) 

leads to: 

^0''^ d^^mp^f + {d,p,f (17) 

-Htanh^ pi{df,ipif + t&nh^ P2(a^'/?2)^ + P2)} 
where: 

P2) = ^ Veff{ps,hs) + 2Mq sinh^ pi sinh^ p2 

S 

(18) 
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This potential can develop one or two minima, depending 
on the values of hg, thus generating either one or two 
gapless modes. 

From the exact solution one can see that massless 
modes appear when h exceeds certain threshold value 
M, where M is a function of the coupling constant 7 
and parameter Mq ( a similar condition is obtained for 
the type II models). In view of the above considerations 
these modes can be interpreted as fluctuations of the su- 
perconducting phase and naively one may identify h = M 
with the onset of superconductivity. 

This interpretation is incorrect however, since it does 
not take into account vortices. Massless phases of the 
Fateev's models, like any other two-dimensional critical 
theories with U(l),or[/(l)x[/(l), symmetry, may be un- 
stable with respect to vortices. The reason being that the 
naive approach doesn't take into account non-analytic 
configurations of the fields, that on a lattice may give a 
finite contribution to the free energy. The real transition 
is of the Beresinskii-Kosterlitz-Thouless (BKT) type ||] 
and occurs at temperature below the mean field transi- 
tion temperature established by the condition h = M. 
One simple way to see the origin of the vortex configura- 
tions of the order parameter field is the following. Trans- 
formations (|l^) , jl^ ) violate an important property of the 
original model (l7|), namely the periodicity of its action. 
The order parameter fields of the original models (||), (||) 
are periodic in ips, while this periodicity is lost after the 
transformations (p^), (|l^). To recover the original peri- 
odicity one should add to the forms (p^,(p^) exponents 
of the dual field (see for instance Q). These terms are 
not contained in the models we consider, therefore the 
latter ones can provide an adequate description of the 
superconducting phase only below the BKT transition 
temperature, where vortices are irrelevant. 

In the next sections we will use the exact solution to 
study the relevance of vortices in our models. To do this 
we shall need to study more carefully the gapless state 
and extract from the Bethe Ansatz equations the scaling 
dimensions of the vortex operators. 



III. THE 

BEREZINSKII-KOSTERLITZ-THOULESS 
TRANSITION IN THE MONOLAYER MODELS 

Let us consider the type-I models first. The Bethe 
ansatz solution deals with the Minkovsky version of the 
models. In our analysis we shall benefit from the fact 
that the ground state energy of the (1-1- l)-dimcnsional 
field theory with coupling constant 7 is equal to the free 
energy of the 2-dimensional classical theory with temper- 



ature T = 7"^ 



F/T^£ 



(19) 



In the Minkovsky version the appearance of the gapless 
state is related to the creation of a condensate of kinks. 
Then using the Bethe ansatz solution one can express the 
ground state energy per unit area in terms of solution of 
the integral equation: 



£(61) - / de'R{e - e')e{e') = m cosh e-h 



F/T = £ = 



M 
2^ 



(20) 



(21) 



where €{9) < for \9\ < B and the integration limit B is 
defined by the condition 



e{±B) = 0. 

The kernel R is related to the two-body ^-matrix: 
1 dlnS'(6i) 



R{e) 



27ri de 



(22) 



(23) 



According to [|l| its Fourier transform, i?(a;), has the fol- 
lowing form: 

1 - R{uj) = - 9)I<A cosh[7rcj(a + 2g)/2h] ^^^^ 



where 



cosh(7ra; /2) smhijTLU / a) 



47r HY + 47rG 

a ■ 



j2 _|_ 4^ 



An + 7^ 



(25) 



with G, H depending on the model. For model (|^) one 
has G — H — 2(n -1-2); G,H always scale with n when 
n ~f 00. It was shown by Fateev that 



27r[l - R{u = 0)] 



g{h/M) 



(26) 



where g{oo) = 1. 

Another quantity that will be useful in the following is 
the dressed charge C(^) jllillll defined by the following 
integral equation: 



C(^) 



d0'R{e~e')c{0') = 1. 



(27) 



-B 



The kernel R is defined in (g4|) and the limit B is deter- 
mined by condition (p^). 

As we have said, two-dimensional U(l)-symmetric crit- 
ical points can be unstable with respect to the pres- 
ence of vortex configurations of the order parameter field. 
The vortices constitute potentially relevant perturbations 
which appear in the effective action as exponents of the 
dual field (see, and, for example , Q). This opera- 
tors become irrelevant if their scaling dimension, c?e; is 
greater than two. When a critical theory includes just 
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one U(l) field a, the scaling dimension of the order pa- 
rameter A = exp(ia), d-A, is related to the scaling di- 
mension of the vortex perturbation as follows: 

dAde = 1/4 (28) 

Hence the superconducting regime exists at 

dA < 1/8. (29) 



The scaling dimension for the primary field $ is defined 



as: 



A+ 



(30) 



where the conformal dimensions determine the 

asimptotics of the correlation functions of primary fields: 



< $A±(2;,0*A±(0,0) >= 



cxp(— 2il?pi?a;) 



{x - wt)2'^+(a; + ivt) 



2 A- 



(31) 



Here 22? is the momentum of the state in units of the 
Fermi momentum p p . 

In (l+l)-dimensional critical theories the scaling di- 
mension of primary fields are related to the finite size 
corrections to the ground state energy For models 

which, as the model in question, has central charge c—1, 
conformal dimensions are given by: 



2A=^(AiV, D, N^) = 2N^ + ± ZD 



(32) 



The quantum number AA'' is characteristic of the local 
field under consideration; in the context of this model it 
represent the number of particles produced by the pri- 
mary field in consideration. The quantum numbers D 
and N^, generate the tower of excited states, and rep- 
resent respectively the number of particles that undergo 
back scattering processes from one Fermi boundary to 
the other and the number of particles added at B{N'^) 
or -B{N-). While AiV and D are fixed by the local 
field, must be chosen to give the leading asimptotics 
in the correlation functions, which is equivalent to mini- 
mize A^ . In integrable models the quantity Z is related 
to the dressed charge introduced above in the following 
way: 



z^as). 



(33) 



For the order parameter operator we have (see Ap- 
pendix) : 



dA = A+(l, 0, 0) + A-(l, 0, 0) = 1/[2Z] 



(34) 



In order to study the stability of our models with re- 
spect to vortices we then have to calculate the value of 
the dressed charge at the Fermi point. We can make two 



preliminary observations. First, &i B ^ oo, correspond- 



ing to h/M ^ oo, we can approximate Eq. ( |27[ ) by the 
Wiener-Hopf (WH) one. Then we have: 

dA = \[l-R{^ = Q)] = \^^ (35) 
4 4 7^ + 47r 

This gives us the upper estimate for existence of the su- 
perconductivity: 



7^ < 47r, 5 > 1/2 



(36) 



This result holds for a^/type I models. Second, since R{6) 
is a non-singular kernel, at small B we have C,{B) — > 1 
and the scaling dimension is too large for the supercon- 
ductivity to occur. Therefore there is a line in the 7 — bi- 
plane separating the superconducting and the disordered 
regions. 

As we noted before for model ( p|) G and H scale like 
n for large n. Then the kernel (|24|) becomes a delta 
function in the limit n — > 00. In this case one has to 
be a bit careful with determining Z^ but the outcome 
is simple: the scaling dimension is /i-independent and is 
equal to 



(37) 



dA = (1 - ff)/4 



We can study in more details the behavior of Eq. (|27 
in the regions _B 3> 1 and B and obtain the asymp- 
totics of the line separating the superconducting from the 
disordered region. The asymptotic form of the scaling 
dimension at large B can be found from the generalized 
Wiener-Hopf method ||l^ and is equal to 

1 



dA = -{l-g)[l + a{M/2hY + ...] 



2h 



(38) 



h + 2g 

Atg— 1/2<<1 the condition dA < 1/8 gives us 

{h/2M)> f{N){g~ 1/2)-^/^ (39) 

with fi approaching 2atA^— >oo. Atn = we have 
fi = 2/il + g). 

To calculate Z for i? ^ 1 it is better to use the follow- 
ing identity: 



2Z = irvp- 



'dh 



where Df is: 



= / dea{e), 

l-B 



(40) 



(41) 



vp is the Fermi velocity and (7(9) is the ground state 
density for e{9). For small B all objects in equation ( |40| ) 
can be calculated explicitly and one has 
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dA — 



1-.9 



M 



M 



(42) 



We have solved Eq.(g^ numerically for various values 
of n and constructed the phase diagram of the model us- 
ing Eq.(p9[) to identify the superconducting region. The 
results are presented in Fig.l. It is clear that the in- 
sulating substrate strongly affects the superconducting 
transition, in particular with increasing n the distance 
between the BKT transition temperature and the mean 
field transition temperature decreases. This means that 
the effective critical temperature increases with n. To 
our knowledge this is the first model that displays such 
characteristic behavior. 
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FIG. 1. Phase diagram of the type I models. The curves 
separate the superconducting region from the disordered one. 



IV. EFFECT OF THE INTERACTION BETWEEN 
SUPERCONDUCTING ORDER PARAMETERS: 
THE TWO-LAYER MODEL 

The situation is considerably more complicated for the 
case of two-layer models. Since we want to focus our 
attention on the effects of the interaction between the 
order parameters we consider only the n = version of 
the model. We have noticed elsewhere that in this case 
the two gapless modes have the same velocity and 
hence at this point the model is conformally invariant 
with the central charge c=2. This result is not general 
for the Fateev's models but valid only for the specific 
case n = 0. However, the results presented below can 
be generalized for a case when the two velocities are dif- 
ferent. In the latter case the low energy sector is split 
into two quasi-independent sub-sectors. The system as a 
whole is not conformally invariant, but each sub-sector 
is (with its own velocity) . Conformal symmetry of each 



sub-sector makes it possible to generalize the basic rela- 
tions between the critical exponents and finite size scaling 
amplitude that we will use in the further discus- 

sion. 

The low-lying excitations for the type-II model with 
rt = are described by the following equations: 



.(+) 



4+j + {1 - i?.(c^)} * = 4. + {^i(^)} * (44) 

where we have used the shorthand notation for convolu- 
tion: 

(• + 00 

{g(c.)}*/-g*/((?)= / dB'g[B -&)!{&) (45) 



{l-i?i(c.)}*6^)=e;(0) 

-f{ifl(c.)}*[4-)+6y] 



(43) 



(-) 



As usual, = ex{d) for 16*1 > B and zero otherwise, 

while = 61(6*) for \B\ < B, and similarly ei+j(A) = 

e±,(A) for |A| > Q and e^^J{X) = e±,(A) for |A| < Q. 
This relations also fix B and Q. The bare energies are: 

e° (0) = M cosh(e') - i (/!+ + /i_ ) , el^^h± (46) 

and the quantities h± are related to the chemical poten- 
tials: 



h± — hi ± h2 



(47) 



Here hi and /i2 must be chosen such that h± are always 
positive. The Fourier transforms of the kernels have the 
following form: 



Ki{n) 



Km ^ - (48) 



sinh(j/f2) ' 



sinh(i^i7) 



^ ' sinhi^^f^sinhf^i. 2A^ ^' 

1-5 

and the remaining constants are: = 1 — g, A = 2 — g. 

For ^{h+ + h-) > M the mode ei becomes gapless and, 
depending on the relative values of h+ and /i_ , it can also 
induce a second gapless mode. In the following we will 
consider /i+ ^ M . In this case e^, is always positive 
and decouples from the other equations such that they 
become: 

ei{e) e^O) + r d9'Ri{e - e')ei{e') (49) 

J-B 

rQ 

+ / dXKiie - \)e_,{X) 
J-Q 

rQ 

e_,(A) = eO_^+ / dAX(A- A')e^.(A') (50) 
J-Q 

i-B 

d9Ki{X - 0)ei{e) 
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For /i_ ^ /i+ it is convenient to invert the kernels in 
Eq.®: 



sinh i^fl 



2sinhrJcosh[(z^ — l)f2] 



.(+) 



+ eiJ^i^h^/2 (51) 



2cosh[(iy- 
Using this form in Eq.(E^) we get: 



* e 



{-) 



.(+) 



X*eP ^ Mcoahe - 
1 



(52) 



.(+) 



where 



2cosh[(iy - 
sinh il cosh[(3i' — 1)17] 



2 sinh M cosh[(i' — cosh[(i^ + 1)11] 



(53) 



Since e^_I^J is very small, to first approximation this re- 
duces to: 

ei(0) = £0(0) + / de'Ri9~d')e,i0') (54) 

where the new kernel R{fl) is given by: 

R{n) = 1 - K{n) (55) 
Clearly become gapless for /i_ < he defined as: 
1 



hr 



dO- 



J-B cosh ^ + cos ■ 



MO) 



(56) 



Already at this stage some interesting features emerge. 
Namely, in order to get the phase transition on the mean 
field level, one does not need both fields /ii, to exceed 
the critical value for a single superconductor. Instead 
of the condition \hi\ ~ |/i2| ~ M/2 we get a weaker 
condition hi > M/2 and h2 > hi — he- 

Now let us discuss how the mean field picture is mod- 
ified by the vortices. For h_ > /ij, when there is only 
one gapless mode we can calculate the scaling dimension 
of the order parameter using the procedure of the previ- 
ous section. However, for /i_ < he when there are two 
gapless modes this procedure should be modified. In this 
case one needs to introduce the dressed charge matrix 



( Zii 


Zi, ^ 




\ Z,i 


Zuv } 


H 



Cii(s) Ci.(Q) 



(57) 



(we follow the convention introduced in ref. ||lj,|l2]) and 
the function ((9) is determined by the following equa- 
tions: 



Cii(^) = 1+ / d9'Ri{e-9')(ii{9') 

J-B 



d\Ki{9-\)Ci,W 



Ci.(A) = / dAX(A- A')Ci.(A') 



d9Ki{X-9)Cii{9) 



= / dAifi(0-A)U(A) 



/ A9Ri{9 - 9')U{9') 

J-B 



U(A) = 1+/ d9Ki{X~9)U{d) 



(58) 



(59) 



(60) 



(61) 



dAi?,(A'- A)C.(A') 

Q 

We note that the structure of this equations is very sim- 
ilar to the one for the Hubbard model away from half 
filling, where the two gapless modes correspond to spin 
and charge excitations [p^ 16|. 

In this general case the formula for the conformal di- 
mensions are given by: 

2A^ (AN, D, N±) = 2Nf + {ZnDi + Z^iD^ 

2 



(62) 



± 



2detZ 



2A±(AN, D, N±) = 2iV± + {Zi.Di + Z,,D, 



(63) 



± 



ZiiAA^^ - Z^iAiVi 
2detZ 



The quantum numbers AN = (AiVi, AiV^), D = 
{Di,D^) and N='= = {Nf,N^) are the obvious gener- 
alization of the one defined in the previous section. For 
two coupled Gaussian models with a total central charge 
c — 2, the scaling dimensions for primary fields are given 
by 



d(AN, D) = [ZiiDi + Z^iD^Y 



{Zi,Di+Z,,D,Y 



Z^^AiVi - Zi^ATV^ 
2dctZ 

ZiiAiV^ - Z^iATVi 
2detZ 



When two sectors of the Gaussian model have different 
velocities vi and V2 , the correlation functions of primary 
fields are given by: 

< $A±(a:,t)$A±(0,0) >= (65) 

exp [-2i{'DiPFi + 'D2Pf2)x] 

(x - ivit)'^^t {x + ivity^i {x - iv2tf^i {x + iv2tf^^ 

The vortices generate operators represented by exponents 
of the dual phases of the two gapless fields. We shall 
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call the scaling dimensions of this operators doi and de„ 
respectively. The conditions for irrelevance of vortices 
are : 



> 2 , de„ > 2 



(66) 



To see the effect of the interaction between the order 
parameters on the stability of the superconducting state 
with respect to vortices we shall consider the dependence 
of the scaling dimension of the potentially relevant oper- 
ators as a function of for fixed 3> M . In particular 
we will compare the scaling dimensions in two extreme 
limits where Eqs (|5^|6ll) can be studied analytically. The 
limit of small which corresponds to the physical sit- 
uation in which the two superconductors have the same 
bare critical temperature, and the limit > he- Having 
fixed /i+ correspond to fix, for example, hi which is pro- 
portional to the mean field critical temperature of layer 
1. We want to observe how the presence of the other 
superconductor affects the scaling dimension of the or- 
der parameter of this superconductor, i.e. his effective 
critical temperature. 

For /i_ ^ 1 the matrix problem given by Eqs. ( |58[]6l| ) 
can be reduced by to scalar one by Fourier transforming 
with respect to A. In these circumstances the dressed 
charge matrix has the form: 



Z = 



f C{B) 



where (^(9) is determined by the following equation: 
({9) = 1+ / de'R{9-9')C{9') 

J-B 



(67) 



(68) 



and the kernel is defined by Eq.(p5|). In this limit the 
scaling dimension is given by: 



d(AN,D) = {aB)Di + ^D,f + '^Dl 
and, as shown in Appendix: 

=d((0,0),(l,0)), 
de„ =d((0,0),(l,-2)). 



(69) 



(70) 
(71) 



Then in general the model will present two different BKT 
temperatures. 

In the limit i? >> 1 Eq.(|68|) can be solved analytically. 
Using again the WH method we obtain: 



as) 



'2v 



from which: 



1 



a 



(72) 



(73) 



Then in this particular limit the conditions ( |66| ) are al- 
ways satisfied and the system is always stable with re- 
spect to vortices. We notice that there are operators 



characterized by Di = 1 and Dy = — 1 or £>i =0 and 
Dy = 1 with a smaller scaling dimension: 



d{{0, 0), (1, -1)) = d{{0, 0), (0, l)) = iy= df 



(74) 



Again following the results in the Appendix one can 
easily see that these operators are associated with the 
bosonic exponents containing linear combinations of the 
dual phase of the two fields. These operators become 
relevant for g < 1/2. 

For h- > he one again has only one gapless mode and 
repeating the procedure of the previous section we get, 
for B > 1: 



de^{h > he) = V, 



(75) 



in agreement with the results obtained for the type I 
models. 

Then the effect noticed at the mean field level survives 
a deeper analysis. The effective critical temperature of 
a superconductor is enhanced by the presence of another 
one. The effect is present only if the critical temperature 
of this superconductor is above some critical value and 
is maximal when the two superconductors have the same 
critical temperature. 



V. DISCUSSION 

We have studied some properties of two families of in- 
tegrable models describing two dimensional superconduc- 
tors interacting with an insulating substrate. The models 
show some interesting features that may be general prop- 
erties of layered superconductors. Due to integrability of 
the models in question one can study them in a strong 
coupling regime. So most of the results presented here 
are non-perturbative, as suggested by the previous anal- 
ysis lul. 

The models presented here can be grouped into two 
families. The members of one family describe a single 
superconducting plane interacting with n insulating lay- 
ers. This gives a possibility to study the dependence 
of the BKT transition temperature on n. Remarkably 
the presence of the insulating stratum makes the system 
more stable with respect to the BKT transition and the 
transition temperature increases with n. The models of 
the second family describe an insulator sandwiched be- 
tween two superconductors, and from this it is possible 
to extract information about the effect of the interac- 
tion between the superconducting order parameters. We 
have considered only the simplest model of this family, 
given by n = 0. Even in this case the presence of the 
other superconductor stabilizes the system with respect 
to the BKT transition. For most of the results presented 
here it is essential that the systems we consider are two 
dimensional. 
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APPENDIX 

In this appendix we summarize some basic results on 
the correlation functions of exponents of bosonic fields 
in Gaussian field theories and present a simple way to 
identify the operators that correspond to various combi- 
nations of quantum numbers AN and D. Let us consider 
the Gaussian model: 



(76) 



It is a well known that the correlation function of expo- 
nents of bosonic fields is given by: 



(exp[i/3i<i>(a)]...exp[i/3jv$(eiv)]) 



(77) 



where z = r + ix/v, z = t — ix/v and in an infinite system 
this is different from zero only if: 



(78) 



The expression (|7^) is factorized into analytic and anti- 
analytic parts and then it is useful to rewrite it as: 



(exp[i/Ji<i>(ei)]...exp[i/3Ar$(^^)]) 

G{zi, ...,Zm)G{zi, ZAr)(5^ p 



where 



c(w)^n(?) 



(ftft74^) 



This factorization guarantees that analytic and anti- 
analytic parts of the correlation functions can be stud- 
ied independently. Since factorization of the correlation 
functions is a general fact, it can be formally written as 
factorization of the corresponding fields. Then inside the 
(...)-sign one can rewrite $(z, z) as a sum of independent 
analytic and anti-analytic fields: 

cf>(z,z) = 0(z) + ,^(z), (80) 
exp[i/3$(z, z)] = exp[i/30(z)] exp[i/3(;()(z)] (81) 

For many purposes it is convenient to introduce the 
'dual' field Q{z, z) defined as 



e(z,z) 



(82) 



(83) 



In order to study correlation functions of the analytic 
and the anti-analytic fields, we define the fields 



A(/3,z)=exp<i-/3[$(z,z)-)-e(z,z)] 
A(^,z)=exp<j^/3[$(z,z)-e(z,z)] 



(84) 



with, generally speaking, different (3, (3. With the opera- 
tors A{(},z),A{P,z) one can expand local functionals of 
mutually nonlocal fields <I> and 6. Suppose that F($, 6) 
is a local functional periodic both in $ and Q with the 
periods Ti and T2, respectively. This functional can be 
expanded in terms of the bosonic exponents: 



F($, e) = ^ F„,™ exp[(2i7rn/Ti)$ + (2i7rm/r2)e] 

71, m 



where 



n 
n 



m 

m 
T'2 



(85) 

(86) 
(87) 



It turns out that the periods Ti , T2 are not arbitrary, but 
related to each other. The reason for this lies in the fact 
that the correlation functions must be uniquely defined 
on the complex plane. We can see how this argument 
works using the pair correlation function as an example: 



(^9) (A(/3„„,zi)i(A 



zi)A{-l3 z2)A{-(3 

nm ; 



1 



\Zl2 



2d 



Z12 
Z12 



where we introduce the quantities 

d = A+ + A- = — 
Stt 



and 



which are called the 'scaling dimension' and the 'confor- 
mal spin', respectively. 

The two branch cut singularities in Eq. ( |8^ ) cancel 
each other and give a uniquely defined function only if 



25* = (integer) 



(89) 



i.e., physical fields with uniquely defined correlation func- 
tions must have integer or half-integer conformal spins. 
This equation suggests the relation 
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47r 



Ti = — EE VAttK 
T2 



(90) 



as the minimal solution. Here we introduce the Luttinger 
liquid parameter K for future convenience. The normal- 
ization is such that at K ^ 1 the periods for the field 
$ and its dual are equal. The quantities A+,A~ are 
called 'conformal dimensions' or 'conformal weights'. In 
the case of the Gaussian model (^6|) the conformal di- 
mensions of the basic operators are given by: 



A+ = 0^ 

nra — A nr: 



?7r = — ( mV K 
8 



1 



n 
K 
n 



(91) 



This equations can be rewritten in terms of the quan- 
tum numbers AN and D previously introduced as: 



2A, 



± 

nra 



Dy/K AN 



(92) 



In integrable models the parameter K is related to the 
dressed charge introduced in Sec. Ill: 



K = 4C(S)2 
from which we obtain the form (|9, 



(93) 



(94) 



Comparing it with equation ( p5| ) we can easily see that 
the scaling dimension 0?$ of the field <& is given by AiV = 
1 and D = 0, while the one of the dual field, rfe, by 
AA^ = and £1 = 1: 



= l/4C(i?)2 



(95) 
(96) 



All the considerations above can be generalized to the 
case of n gaussian fields with the same velocity: 



S=-J drdx^[i>-i(9,$„)2 + i;(9,$„)'] 



(97) 



In a similar fashion as before we can introduce the fields: 



A(B,; 



exp<j -B-[*(z,z) + 0(z,z)] 



(98) 



^(B, z) EE exp <j -B • [*(z, z) - 0(z, z)] 



where * =_ ($1, $„), ® = ei,...,e„, B = (/3i, /?„) 
and B = {(3i, ...,/?„). The generalization of the conformal 
dimension and spin are: 



d(AN,D) = l/4AN^X-iAN + D^XD = A+ + A" 
S'(AN,D) = AN^D 

where AN = (AA^i, AiV„)^, D = {Di,...,Dnf and 
X = Z with Z is the dressed charge matrix defined 
by (|5^). Then the conformal dimension of each field is 
given by: 



A+ = 



1 



(99) 



For the case of two gapless modes of interest in the paper 
the forms (^ and ( p9| ) reduces to 

A(B,z)A(B,z) = (100) 
exp [i{ZiiDi + ^21152)61 + i{Z22ANi - Zi2AN2)^i 
+ i{Zi2Di + ^21^2)62 + i(^iiAiV2 - Z2iAiVi)$2] 



It is easy to read from ( |100D that the exponent of the 
dual field 81 is characterized by Z12D1 + Z21D2 = and 
AiVi,2 = 0, while for 62 you need ZnDi + Z21D2 = 
and ANi^2 = 0. Then for a dressed charge matrix of the 
form (|6^) the conformal dimensions for exponents of the 
two dual fields are 



dei =d((0,0),(l,0)) 
de, =d((0,0),(-l,2)). 



(101) 
(102) 
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